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Abstract: In this study, nonlinear singularly perturbed problems with nonlocal condition are evaluated by finite 
difference method. The exact solution u(x) has boundary layers at x = 0 and x = 1. We present some 
properties of the exact solution of the multi-point boundary value problem (1 )-(3). According to the perturbation 
parameter, by the method of integral identities with the use exponential basis functions and interpolating 
quadrature rules with the weight and remainder terms in integral form uniformly convergent finite difference 
scheme on Bakhvalov mesh is established. The error analysis for the difference scheme is performed. £ — 
uniform convergence for approximate solution in the discrete maximum norm is provided, which is the first- 
order (O(h)). This theoretical process is applied on the sample. By Thomas Algorithm, it has been shown to be 
consistent with the theoretical results of numerical results. The results were embodied in table and graphs. The 
relationship between the approximate solution with the exact solution are obtained by Maple 10 computer 
program. 

Keywords: Singular perturbation, finite difference method, Bakhvalov mesh, nonlocal boundary condition, 
uniform convergence. 

1. Introduction 


In this paper we shall consider nonlocal singularly perturbed nonlinear problem 

— £ 2 u" + fix, u) — 0 , 0 < X < 1, 
u( 0) = 0 


i 

k 0 u( 1) = / fejUfSj) + k m+1 I u{x)dx + d 


(1) 

( 2 ) 

(3) 


f{x,if) e C 1 ([0,l]xf?), 


df 

du 


ix,u) > a > 0, Sj £ (0,1), i = 1 , k 0 > 0 


where, 0 < £ « 1 is small perturbation parameter. 


Singularly perturbed differential equations are symbolized by the existence of a small parameter £ — multiplying 
the highest order derivatives [4], Like these problems, thin transition layers where the solution varies rapidly, 
while away from layers it behaves regularly and varies slowly. Classical numerical methods are inappropriate 
for the singularly perturbed. So, it is important to develop suitable numerical methods to these problems, whose 
accuracy does not depend on the parameter value £; that is, methods that are convergence £ uniformly [18]. 
Such problems arise in many areas of applied mathematics: fluid mechanics, chemical-reactor theory, control 
theory, electrical networks, and other physical models. In recent years, Singularly perturbed differential 
equations were studied by many authors in various fields of applied mathematics and engineering. For 
examples, Cziegis [1] studied the numerical solution of singularly perturbed nonlocal problem. Cziegis [2] 
analyzed the difference schemes for problems with nonlocal conditions. Bakhvalov [3] investigated on 
optimization of methods for solving boundary-value problems in the presence of a boundary layer. Amiraliyev 
and Cakir [5] studied numerical solution of a singularly perturbed three -point boundary value problem. 
Amiraliyev and (,'akir [6] researched numerical solution of the singularly perturbed problem with nonlocal 
boundary condition. Amiraliyev and Duru [7] studied a note on a parameterized singular perturbation problem. 
Boglaev [8] applied uniform numerical methods on arbitrary meshes for singularly perturbed problems with 
discontinuous data. Adzic and Ovcin [9] obtained nonlinear spp with nonlocal boundary conditions and spectral 
approximation. Amiraliyev, Amiraliyeva and Kudu [13] analyzed a numerical treatment for singularly perturbed 
differential equations with integral boundary condition. Herceg [10] estimated on the numerical solution of a 
singularly perturbed nonlocal problem. Herceg [11] analyzed solving a nonlocal singularly perturbed problem 
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by splines in tension, fakir [14] researched uniform second -order difference method for a singularly perturbed 
three-point boundary value problem. Geng [12] gave a numerical algorithm for nonlinear multi-point boundary 
value problems, fakir and Arslan [15] applied a numerical method for nonlinear singularly perturbed multi¬ 
point boundary value problem. Bitsadze and Samarskii [16] studied on some simpler generalization of linear 
elliptic boundary value problems, fakir [17] obtain a numerical study on the difference solution of singularly 
perturbed semilinear problem with integral boundary condition.fakir [5] applied the difference method on a 
Shishkin mesh to the singularly perturbed three-point boundary value problem. 

Some methods for singularly perturbed problems have been studied in different ways [9-12]. In this 
study we present uniformly convergent difference scheme on an equidistant mesh for the numerical solution of 
the problem (1)—(3). The difference scheme is constructed by the method integral identities with the use 
exponential basis functions and interpolating quadrature rules with the weight and remainder terms in integral 
form [18]. In Section 2, we obtain some analytical results of the multi-point bounday value problem (l)-(3). The 
difference scheme constructed on Bakhvalov mesh for numerical solution (1)—(3) is described in Section 3. We 
prove that the method is first-order convergent in the discrete maximum norm. In the following Section, a 
numerical example is considered and constructed which is in agreement with the theoretical results. According 
to results, the finite difference method on Bakhvalov mesh is more powerful method than other methods for 
nonlinear singularly perturbed multi-point boundary value problem. 


2. The continuous problem 

Here we give some properties of the solution of ( l)-(3) with Lemma 2.1. Because, they are needed in 
later sections for the analysis appropriate numerical solution. We use Halloo for the continuous maximum norm 
on the [0,1], where g(x), is any continuous function. 

Lemma 2.1. 

Let f{x, u ) £ C 1 (\0,l]xR) and — (x, u) is uniformly bounded in u(x). We assume that 


m 

^/qwoOi) 

i=1 



< K 


(4) 


where w 0 (x) > |w(x)|, 


-e 2 w' + a(x)w(x) = 0 

(5) 

c> 

II 

p 

II 

(6) 


the solution of this problem. 

So, the solution of the Eq. (1) satisfies the inequalities 

IK*)llc[o,i] ^ C 0 (7) 


and 


\u (x)| < C jl + - 


—sfax —sfail—xy 
q e G £ 


•,0 < x < 1 


where, C 0 and C are constants independent of £ and the h. 


( 8 ) 


Proof. We use intermediate value theorem for /(x, u) in the Eq.(l) and we obtain the following the equalities, 


fix, u) — fix, 0) df 

- — ——ix,v), v = vu, 0<y<l 

u - 0 ou 


www.ijmret.or 


g 


Page 26 










( \ 

International Journal of Modern Research in Engineering and Technology (IJMRET) 

www.ijmret.org Volume 1 Issue 5 II December 2016. 

, s s df 

fix, u ) = f{x, 0 ) + u(x) — (x, i9) 
ou 


= —Fix) + a(x)u(x) 


where a(x),F(x) are sufficiently smooth on [0,1] and \Q,l]xR. 


Also, 


df 

— (x, i9) = a(x) > a > 0 , i9 — intermediate value, 

ou 


So, we obtain the following linear equation, 


—e 2 u" + a(x)u(x) = F(x), 0 < x < 1. 

(9) 

Now, we can write the solution of the Eq.(9) as follows 


u(x) = vix) + 2 w(x) 

( 10 ) 

—e 2 v'' + a(x) v(x) = F(x) 

( 11 ) 

o 

II 

o 

II 

o 

( 12 ) 

where v(x) is solution of the Eq. (11)- the Eq. (12). 


First, we give the estimate vix), 


b(x)| < 117(0)1 + |t 7 (l)| + a _ 1 ||F|| tx) < C 1 . 

(13) 

Second, we give the estimate w(x). 


|w(x)| < |w( 0 )| + |w(l)| + cc —1 HOlloo < 1 

(14) 

Then, from the Eq.(13) - the Eq.(14) we have the following inequality 


|u(x)| < b(x)| + Ul|w(x)| < C x + m 


|u(x)| < C 0 . 

(15) 

we now prove the estimate the Eq. ( 8 ). 


if u (x) is pulled from the Eq. (9), we obtain 


-1 

u (x) = — 5 - [F(x) — a(x)u(x)] 

£ z 

(16) 

and from the Eq. (16) 
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\u" (x)| < 


■ [F(x) — a(x)u(x)] 


<—■■ 
£ z 


Now, we take derivative of the Eq. (9) and if it called u (x) = v 0 (x), 
the Eq. (9) takes the form with boundary condition 

-e 2 v 0 " + a(x)v 0 (x) = 0(x) 
u'(0) = v 0 (0) = 0,u (1) = v 0 (l) = 1 

In a similar manner, we proceed with the estimation of 0(x), u (0) and u (1) respectively, 
from the Eq. (7), we obtain the following equation 

0(x) = F (x) — a (x)u(x) < C 1 . 

We use the relation lor g 6 C 2 as follows: 


9 00 = 


g(P)-g(a) 


ft — a 


s 

J a 


P-t 


ft — a 


-T 0 (x- t) 


g" (t'tdt, a < x < ft, a < ft 


where 


r 0 (x-t) = [ 


1 , x — t > 0 

0 , x — t < 0. 


the Eq. (21) with the values g{x) — it(x), a = 0, /? = £, x = 0 and from the Eq. (7)- the Eq. (17) 


|u'(0)| < 


u(e) — u( 0) 


+ 


f |u"(t)|dt 
Jn 


< 


(17) 


(18) 

(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Similarly, 

the Eq. (21) with the values g(x ) = u(x), a = 1 — s, ft = 1, x = 1 and from the Eq. (7)- the Eq. (17) 


|u'(l)l < 


u(l) — u( 1 — £) 


+ 


f 1 I u" 
Jl—£ 


(t)\dt < 

£ 


(23) 


We write the solution of the Eq. (18)- the Eq.(19) in the form, 

»oW = bW + v 2 (x) 

where v 1 (x), v 2 (x) are respectively the solution of the following problems, 

—e 2 v 1 " (x) + a(x)v 1 (x) = 0(x) 
tb(0) = 0,u'(l) = Vi(l) = 0 


(24) 

(25) 
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—£ 2 v 2 " (x) + a(x)u 2 (x) = 0 , (26) 

V 2 (0) = bo(0)|,u'(l) = v 2 (l) = bo(l)| 

According to the maximum principle in the Eq.(24)- the Eq. (25),we can give the following Barrier function, 

t/'iCx) = +vfx) + a _1 ||0(x)IL 

This Barrier function provides the conditions of the maximum principle and 

v 1 (x) < C. 

In a similar manner, according to the maximum principle in the Eq.(26), we can write 

u 2 (x) < 0 (x) 

where 0 (x) is the solution of the following problem with constant coefficient, 

—e 2 6 (x) + a 0 (x) = 0 


0(0) = K(0)| < - , 0(1) = K(1)| < - 

£ £ 


(27) 

(28) 


where a(x) > a > 0 and the solution of 0 (x) as follows, 


0(x) = 


V n (°) 

e £ — 

fa{l-x)' 

e £ 

+ v„(l) 

—fax fax' 

e £ — e £ 


•fa —fa' 

—e £ + e £ 



(29) 


after some arragement, we can obtain, 


C 

e{x) < - 

£ 


•fax —yfa(l—x) 


e e + e e 


(30) 


Finally, from u (x) = v (] (x), |u 1 (x)| < C, 0(x) < v 2 (x), v () (x) = u, (x) + v 2 (x), we have the following 
inequality. 


u'(x)\ = b 0 (x)| < \vfx)\ + |u 2 (x)| 

fax —fafl-x)' 


c 

<c + - 

£ 


6 £ 6 £ 


£C ( 1 + I 


•fax —fa(l-x) 


e £ + e £ 


which leads to the Eq. ( 8 ). 


3. Costruction of the Difference Scheme and Non-Unifom Mesh 

In what follows, we give by the following non-uniform mesh on [0,1]: 


( S i- N ) 

= [0 = x 0 < x 1 < ••• < x w _! <x N = l,x N . = Si,N t = —j—j 




We present some properties of the mesh function g (x) defined on a) N , which is needed in this section for 
analysis of the numerical solution. 
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9x.i 


9x,i 


9i 9i—l 


9i +1 9i 




+i 


9x,i 


9i+l 9i 


9xx.i 


hi ' ax ‘ l hi 

9x,i — 9x,i 1 (9i+l ~ 9i 9i ~ 9i—1 


_ 1 ( 9i +1 ' 
hi V 


hi = — X; , hj = 


hi “i v ft i+ i 

K + h i+1 


K 


\\g\L = \\g\L, ^ = max| 5i | 

The difference scheme, we shall construct as following from the identity 

r x i +i :: r x i +i _ 

hj - I — £ 2 u" (x)<pj(x)dx + hj -1 I f{x, u)(p t (x)dx = 0, i = l,N — 1 


'* 1-1 




where {<"/),■ (x)}, i = 1, IV — 1 are the basis functions and having the form 


<Pi (1) (x) = ■ 


x — X, 


i—1 




<Pi0c) = < „(2)<V> -^i±l 


<Pi(x) 

o, 


— X 


h. 


■i +1 


X;^ < X < X; , 

X; < X < X i+1 , 
X g (Xi.^Xf+i) 


<Pi ^ (x) and 1 (x) respectively, are the solutions of the following problems, 

—£ 2 q>" — 0 

<p(Xi_i) = 0, <p(x ; ) = 1 


— £ 2 (p" — 0 


cp(Xi) = l,<p(x i+1 ) = 0. 

If we rearrange the Eq. (31) it gives 

-1 f Xi + 1 n . f Xi + 1 

ij I u (x)(pj(x)dx + h; I f{x,u)(p i {x)dx = Q, 


-e z h ; 


(31) 


(32) 

(33) 

(34) 

(35) 


t = 1,1V — 1 (36) 


i 




and after a sample calculation 


where 


-l f * i+1 - 

£ 2 h; I U (x)<P;(x)dX + /(Xj,Uj) + /?; = 0, 


i = 1,1V- 1 


r*i+i r*i+i ^ 

Jti=hiJ dxtpi(x) I —/(x, 

x i—l x i -1 


u)KZ(x,OdS, i — 1, N — 1 


(37) 


(38) 
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and 


KSb, 0 = t q (x-0-T 0 ( i+1 2 i_1 -f) + (Xt + 1 -- 0 ° ( i+1 2 i-1 -*)■ 

So, from the Eq. (37), we propose the following difference scheme for approximation the Eq. (l)-the Eq. (3): 


—e 2 % £i + /(Xj,Uj) + /?, = 0 t = 1, IV — 1 


(39) 


It is now necessary to define an approximation for the second boundary condition of the Eq. (1). We accepted 
that x w .is the mesh point nearest to s ; . 


m 1 

gt In — / /qt ijy. T k m+ 1 I u{x)dx + d 

ti 1 J ° 


N 


= Y k < UN i + km + i Y hiUi+Ti 

i =1 U=1 


+ d 


(40) 


where remainder term 


N 

V 1 f Xi d 

r i =YJ ^ ~ Xi ~ 1 - 1 dx u ^ d S- 


(41) 

- , ii y 

— x ‘-i 

By neglecting R t and r t in the Eq. (39) and the Eq. (40), we may propose the following difference scheme for 
approximating the Eq. (1)- the Eq. (3): 


- £ 2 y x -x,i + f(x t ,y t ) = 0, i = 1,N 
y 0 = 0 

m N 

= Y kiyN i + k m+l Y hiyi + d 


(42) 

(43) 

(44) 


i=l 


1=1 


We will use the Bakhvalov mesh to be £- uniform convergent of the difference scheme the Eq. (42)- the Eq. 
(44). Transition points of the Bakhvalov mes are defined as 


a = min 




and Xj mesh points as follows: 


a < — , X; G [ 0 , a] , X; = 
<7 = i , Xj £ [0,<x], x,- 
x e [a, 1 — a], Xj = a + 


—{fa) £ln [l — (1 — e) , i — 0,1, 





-(W) 1 
e 


4 i 
N 


,i = 0,1, 


(,_3*cO iI = ! +1 ._ 


3JV ft(1 ) _ 2(1-20-) 
4 ' JV 


W 

4 
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Ci 6 [1 — a, 1 ], x t = 1 — a—{fa) £ln ( 1 — (1 — £) 


1 — a = —, x t e [1 — a, 1 ], Xj = — (fa) eln 



3N . . 

——— , l = 

— + 1 , 

N ) 

4 

-(W) _1 ) 

4-^)1 

1 — e 4 e 


) 

N 


3 N 

i = — +1.....N 
4 


where, N is even number, er « 1. 


4. Uniform Error Estimates 

Let z — y — u, x 6 c5 N , which is the error function of the difference scheme the Eq. (42)- the Eq. (44) and the 
solution of the discrete problem 

Iz = -E 2 z M l + [/(X;,y;) - fiXi.uJ] = R t , i = 1, N - 1 (45) 


z 0 = 0 


feo% = 'Yj k i ZN i + fem +l Y hiZ ‘ =ri ' i = 1 ' 


Nj T n -m +1 
i=1 i =1 


N 


(46) 

(47) 


where R t and r t are defined in the Eq.(38) and the Eq.(41). 

Lemma 4. 2. The solution z, of problem the Eq. (45)- the Eq. (47 ) satisfies following 


\\R\L# h < CN- 1 


< civ- 1 . 


Proof. We shall evaulate reminder term R h 
1) For a < —, X; 6 [0, er] 

4 

For X; = — /? _1 £ In [l — (1 — e) , i = 0,1,..., 
we obtain 


hj — X; — Xj_ t = — (fa) £ln 


4t' 


+ (V«) eln 


1 - (1 - £) 


4(t - 1) 


N 


According to t, we use the intermediate value theorem as follows: 




(48) 

(49) 


and we have 


h L + h i+1 

hj = 1 - i+1 < CN- 1 . 
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-fax 


-V5s i+1 1 £ln [l-Cl-e] 4 ^ 1J ] ) -V^-(W) 1 £ln[l-(l-£> 


4(i+l)T 


e ^ — e e < e 


— e 


< e ln[l-d-,)^] _ e >n[l—(I - < e ln[l-d-,)^] _ 


< e 


ln [_ (1 _ £ )lCi=12]_ ln [_ (1 _ £) iC^]< , 


and 


—V«( l-x i+1 ) -Va~(l-r,-_ 1 ) 

e £ — e £ < CIV -1 


We use above inequalities and we can obtain 

R t < 2Ch t < CN- 1 . 


2) er = —, e [0, er] 

4 

We use the followig equation, 

h t = X; — x i _ 1 = — /? -1 £ln 

we obtain. 


/ -M 4tl 



( zl\ 4 i 

1 - l- e4 £ — 

+/? 2 £ln 

1 - 

1 — e 4 e — 

V J N. 


V / IV J 


r-l 


and 


4(1 — e 4 e )1V 

ft; = /? -1 £--- — -< 4/? -1 lV -1 < CJV -1 

1 - ^4(1 - e^)JV -1 


h, + h i+1 

h; = 1 - i+1 < CN- 1 . 


After doing some calculation, 

R t < CN - 1 . 

3) X; G [er, 1 — a], X; = er + (i — -) h^\ i = —+ 1 ,= 2(1 2a \ x L G [er, 1 — er], 

v 4 / 4 4 N 

Let er = —Basins < - , 

According to t, we use the intermediate value theorem and we obtain, 

hi < CN- 1 

and 

^ = h (1 ) < CIV -1 . 

Consequently, 

i?i < CIV -1 . 

1 

4) x ; G [cr, 1 — cr], er = - 


and 


We have 


hi < CAT 1 
h, = h« < CIV -1 . 
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R. < CN 


-l 


5) xi £ [1 — a, 1], X; = 1 — a + x e In 1 — (1 — e) 


<‘-tT 


N 


3 N 

i = —+ 1 
4 


where, according to i, we use the intermediate value theorem as follows: 


h, = X; — Xj_ x = 1 — er + /? 1 eln 1 — (1 — s) 


Thus, we obtain 


and 


— 1 — a + /? 2 eln 1 — (1 — e) 


h; + h i+1 

h; = ' - i+1 < CN- 1 


*tf)' 

N 

«(<-?)' 


N 


R t < CN- 1 . 


6 ) 1 — er = —, Xj e [1 — er, 1], x, = —/? 1 £ In 
4 


we use the intermediate value theorem and we obtain. 


. (. 3N\ 


N 


3 JV 

,i = — + 1,...,N 
4 


h, + h i+1 

h, = 1 _ 1+1 < CN- 1 . 


R 


-l n +i -i r +i 

i= /Ijhj I Chi<Pi(x)dx + /Ijhi I Ch i cp i (x')dx < CN 

dY: ‘'Xi -1 


^i -1 


After all these calculations, we have 


IlKlUs* < CAT 1 . 


W /4 


r,l< 2 / ' ( 1 + 7 

• - r< " i=i x i-l 


i=l ^-1 


—yfax —sfa{l—xy 


e £ + e £ 


dx 


+ v h^r (i+~ 

i= N U+i JCi - 1 £ 

+ y h®p(i+- 


—yfcix —VcF(l—x) 


e £ + e £ 


dx 


-sfax —sfa(l—xy 


e £ + e £ 


dx 


= 3W / 4 +i 

After a sample calculation for r ; , we obtain, 

r, < CAT 1 . 

Lemma 4. 3. Let z ; be solution of the Eq. (45)- the Eq. (47). Then there is the following inequality: 

NU* < C[||«IL^ + |r|] (50) 

Ispat. If we rearrange the Eq. (45), it gives 
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where 


Iz — —£ i Zxx i i + a t z t i = 1, N — 1 

dF 

a t = — if, u t + yzf) , 0 <y < 1 

ou 


we can use the maximum principle, and so it is easy to obtain, 

\\z\\ m ,w N <\z N \ + a -1 (Halloo.w/v + H) 


< 


m tv 

+fc m+ i2^Mi + d-r t 

i =1 i=1 


+ a 1 (|| J R|| 00 , W j V - + \r\) 


< oc 1 (Il-Rlloo, w N + M) 

< C[||«L |W# + |r|] 


(51) 


Theorem 4.1. Let u(x) be the solution of (1 )-(3) and let y be the solution of Eq. (45)- the Eq. (47). Then under 
the same assumptions of Lemma (4. 2) and Lemma (4.3) the following s —uniform error estimate 


Ily-wlloo.wA, ^ CN llnN 


(52) 


holds. 

5) Algorithm and the Numerical Example 

We apply the scheme the Eq. (42)- the Eq. (44) to calculate the solution of approximation of the 
following problem: 


—£ 2 u" + u 2 (x) — f{x ) = 0, 
u(0) = 0 

u(l) = u(0.5) + d 


(53) 

(54) 

(55) 


which has the exact solution, 


/ iz£ *\ 

(2x — 1) (e £ — e ?) 

u(x) =- 1 -- + 1. 

e7 — 1 

We solve the difference scheme the Eq. (42)- the Eq. (44) using the following iteration technique, 
2 e 2 df 


' £ 2 ' 

v (n) 

6 A. 

yi-i 


Mm dy 


+ ^rX x oyt) 


(n) i 
Vi + 


r> A 


+i 


y-:lm 


= fi Xi ) +/(x j ,y j (n - 1) )= 1. N — 1, n= 1.2,... 


y 0 


(n) 


= 0, ='Y j k i yf i 1] + km+^h^ 11 1] + d (56) 

i=l i =1 
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‘ £ 2 ■ 

v (n) 

Mi. 

yI —1 


2 e 2 df 


KK +1 dy 




(ll) I 

Vi + 


M 


+i 


(n) 

y i+ i 


/M)+ /(x i ,y i (n 15 ) -y[ n 1) ), i = 1 , - - 1 , 

y^ n) = Mn-i. yw n) = Mn-i - 1, Mo = C 0 > 1 


(57) 


where 


(£f y|_ lin + (iS ^.n - 1 + (^l.n-l) 

_ / £ \ 2 

2 fc) + ^n-> 

The system of the (57) is solved by the following procedure, 

2 £ 2 5 / 


A,- = 


B, = 


C, = 


+£(w ) 


1 hff 1 h 1 Mm dy 

r 1 .f f (? 2 —1)A , (n—1) df f (n—1)^ 

F i = -f[. x i>yi J + Vi Qj {xt.yf ) 

a 1 = 0 , ft = 0 

aw = 0, (In = Mn-i 

2 2 

Bt n F i +Af i 

a, + i = —--—, p , +1 = —--—, i = l, ...,1V —1 

+ C t — A^a^ + C i —A i a i 

y [ (n) = « i+ iy^ + Pi+i, y / 05 = o .5 , i = n - 1,... ,2,1. 

It is easy to verify that 

21; > 0, B t > 0 , C t > A t + B it i = 1,... ,1V. 

For this reason, the described procedure above is stable. Also, the Eq. (42)- the Eq. (44) has only one solution. 
In the computations in this section, we will take d = — 1. 

The initial guess in the iteration procedure is yf* = 0. 5. The stopping criterion is taken as 


max 


|y(»+i)_y(»)| <! 0 —5. 


The error estimates are denoted by 


and 


0 N _ 


lly-wlloo, 


,a> N 
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e N — maxef. 

£ 

The corresponding e —uniform convergence rates are computed using the formula 



The 

numerical 

results obtained 

from 

the 

problem of the 

difference scheme 

by 

comparison, 

the error and 

uniform 

rates 

of convergence 

were obtained and 


these are presentedin Table 1. 


Table 1. The computed maximum pointwise errors e N and e 2N , the numerical rate of convergence P N on the 
Bakhvalov mesh co N for different values of N and £. 


e/N 

8 

16 

32 

64 

128 

256 

512 

1024 

2 -io 

0.1844643708 

0.1031176337 

0.0525233947 

0.0264140654 

0.0131517954 

0.0064588555 

0.0031701248 

0.0012507942 


p=0.826 

p=0.973 

p=0.99l 

p= 1.006 

p= 1.025 

p= 1.026 

p= 1.341 


2 -12 

0.2186063415 

0.1221878445 

0.0639603910 

0.0335687579 

0.0176559863 

0.0092080361 

0.0047118041 

0.0023568060 


p=0.839 

p=0.933 

p=0.930 

p=0.926 

p=0.939 

p=0.966 

p=0.999 


2 -14 

0.2468061329 

0.1374469602 

0.0731044491 

0.0393749598 

0.0214435825 

0.0116955952 

0.0063238248 

0.0033579638 


p=0.844 

p=0.910 

p=0.892 

p=0.876 

p=0.874 

p=0.887 

p=0.913 


2-1(5 

0.2701179896 

0.1497045164 

0.0803967387 

0.0439999693 

0.0244759852 

0.0137180219 

0.0063238249 

0.0032650757 


p=0.851 

p=0.896 

p=0.869 

p=0.846 

p=0.835 

p= 1.117 

p=0.953 


2-18 

0.2895653366 

00.1596652402 

0.0862725672 

0.0477070838 

0.0268973496 

0.0153331458 

0.0087696300 

0.0050041978 


p=0.858 

p=0.888 

p=0.854 

p=0.826 

p=0.810 

p=0.806 

p=0.809 


2-20 

0.3059909186 

0.1679963817 

0.0911603753 

0.0506828286 

0.0288403106 

0.0166239446 

0.0096388095 

0.0055920974 


p=0.865 

p=0.881 

p=0.846 

p=0.813 

p=0.794 

p=0.786 

p=0.785 



The results show that the convergence rate of the considered scheme is essentially in accord with theoretical 
analysis. 



X 


Figure 1. Comparison of Approximate solution and Exact solution for £ = 2 14 
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Figure 2. Error distribution for £ = 2 14 ,2 16 ,2 18 using N=256 

From the graps it is show that the error is maximum near the boundary layer and it is almost zero in 
outer region in the Figure 2. Approximate solution compared with exact solution in Figure 1. 
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